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Section A 3 


Answer no more than THREE questions from this section. All questions in 
this section have equal marks. 


Question 1 


A good approximation to the equation governing the voltage u(x, t) along 
a submarine cable is given by 

ao _ , ae 

at ax?” 
where & is a positive constant and x, t represent distance along the cable 
and time respectively. 


Suppose that the sending end of the cable (x = 0) is maintained at a con- 
stant positive voltage vo and that the receiving end (x = /) is earthed (i.e. 
maintained at zero voltage). Show that the steady-state voltage along the 
cable is vọ(1 — x/l). 

After this steady state has been established the sending end is earthed. This 
takes place at time £ = 0. Determine the voltage u(x, t) at subsequent times. 


Question 2 


The explicit finite-difference scheme 


l 1 
z 4; = 7 ny; 


is used in approximating the problem 


2 
wzy 0<x<l,t>0, 


U(x, 0 =x) O<xel, 
UO) =e) 1>0, 
UD) =g) ¢>0. 


If the mesh ratio r = k/h? satisfies the inequality r < 4, show that the 
scheme's solution satisfies the maximum principle 


max uy, < Maxuy,, , 
S Ca 
where, for any positive valuc of T, D, is the set of mesh points interior to 
the solution domain with £ < T and C, is the set of mesh points on the 
initial line and boundaries of the solution domain with ¢ < T. Deduce from 
this that the scheme is stable for r < }. 


What is the main practical disadvantage of this explicit method, and what 
method would you suggest might be preferable in practice? 
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Section A — continued 


Question 3 


Consider the problem 


>- Vu+ku=F inD, 
u=f onc, 


where D is a plane region bounded by the closed curve C, k is a nonzero 
constant, and F, f are functions defined on D, C respectively. 


Use Green's theorem to show that if k < 0 any solution to the problem is 
unique. 


Suppose that 

D= {(x%, y):0 <x <a, 0< y< b} y 
and that 

k w P 

A a 


for some pair of nonzero integers m, n. By finding a solution to the corres- 
ponding homogencous problem, or otherwise, show that in this case the 
original problem cannot have a unique solution. 


Question 4 


Find the Green's function for the system 
u’+us -f O<x<x, 
u'(0) — u(0) = u(x) = 0. 


Hence write down an expression for the solution to the problem and verify 
that this satisfies the given boundary conditions. > 
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Section B m4 
Answer no more than TWO questions from this section. All questions in this 


section have equal marks. 
> 


Question 5 


Verify that the operator 


a a? 
La zat 2xaa t+ T 4) oa +22 + Gx +2) 2 
is hyperbolic. Show that its characteristic coordinates may be defined by 
E=- De, 


n=(x+2)e', 


and hence derive the standard (canonical) form of the equation L[u] = 0° 
as 


Ou | du 
g- n) an + 3 


Find the general solution of this last equation in the form 


ul, n) = SG) + a(n) + (1 — O'E) 


where f, g are suitably differentiable but otherwise arbitrary functions. 


=0. 


Deduce that the solution of the problem 
Liuj = 0 xeR,t>0, 
u(x, 0) = p) | xeR, 


Hx =@-xP'@) xeR, 


ux, 1) = etna) + 4 ple)de , : 


where 
wa = xe? + Al —e7'). 
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Section B — continued 


Question 6 


Show that the finite-difference scheme 


1 1 
k AAC + 052)u,)} = j ny 


approximating the problem 
ôU aU 


7 ae O<x<I,t>0, 
U(0, t) = U(I1, 1) = 0 t>0, 
U(x, 0) = f(x) O<x<l, 

is unconditionally unstable for « > 4sec?4nh, and stable if 
2 +r<sł}, 


wherer = k/h? is the mesh ratio. For what valuc of ais the scheme explicit ? 
Show that the value @ = — 4r gives the Crank-Nicolson scheme for the 
problem. What can you deduce about the stability of cach of these schemes? 


Use the Taylor series expansions 


h? @ ht as 
62U (x,t) = 2 (Fm + a =) U(x, 1) + OA) , 


3 ka 
A,P (x,t) = («3 + Ta) POs t) + Ok’) , 


to show that when a is independent of r the finite-difference scheme has a 
local truncation error of order O(k) + O(h?). Show that the Crank- 
Nicolson scheme has a local truncation error of order Ok?) + O(h?). 
Find also a value of a for which the scheme has a local truncation error of 
order O(k?) + O(kh?) + O(h*). > 
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Section B — continued 


Question 7 


A square ice rink of side / is bounded by a pipe maintained at a steady 
temperature u,°C, where u, < 0. The rink is insulated beneath, but gains 
heaf from the atmosphere above it which remains at a constant temperature 
uC, where u, > 0. 


If the ice is sufficiently thin, the rink may be adequately represented by the 
region [0,/] x [0,/] of R?. By considering the heat balance in a small 
square of the rink bounded by the lines x = Xo, X = Xo + Ax, y = yo, 
Y = Yo + Ay, or otherwise, ‘show that the steady-state temperature 
u(x, y) of the rink may be described by the equations 
2 2 
Fa t ot + My, — u) =0 O0<x<lO0<y<il, 


u(x, 0) = u(x, 1) = u 0 E 

uO, y) = u(l, y) = My 0 l, 

where A is a real constant. You may assume 

(a) that the rink gains heat from the atmosphere at a rate per unit area 
Proportional to the temperature difference between the rink and the 
air above it; 


(b) that the rate of flow of heat per unit arca across any surface within 
the ice is proportional and of opposite sign to the temperature 
derivative normal to that surface. 


Using the finite Fouricr transform method, show that the temperature 
distribution of the rink can be expressed as 


o ý oA 
4P Sa cosh p, (7 5 +) -mmx 
u(x, y) = u + (tia — ty) mpi l- — sin > 
mode 


* cosh tp, 


where p, = n?n? + PA Derive a condition independent of x and y 
sufficient to ensure that the entire rink is in fact frozen. 


` 
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Section B — continued 


Question 8 


Let v(x) be any solution of the equation 
v +Ipv=0 x>0, 


where p is a monotonically increasing positive function of x and 2 is a 
Positive constant. Let a, be the nth non-negative zero of v. By defining new 
coordinates 


Xex-a i=an-l,n, 
and corresponding functions 
px) = p(x), 7 
v(x) = v(x) , 
use the Monotonicity Theorem to show that 
Onti — Ay < An — a, 2 =2,3,..., 
Let 2, be the second eigenvalue of the system 
u" + AL + xu =0 O<x<1, 
u(0) = u(1) =0. x 


By comparing the system with others having constant cocfficients, show 
that 


2R? <2, <4n?. 


It is known from the Oscillation Theorem that the second eigenfunction of 
the system above has precisely one zero (say a) in the interval (0, 1). By 
considering the problem 


u” + XL + xu = 0 O<x<a, 
u(0) = ua) =0, 


prove that 


t<a<VJh. 
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Section B — continued 


Question 9 


Write down the five-point finite-difference approximation for the problem 
W2U(x, y) = f(x, y) O<x<10<ye<l, 
U(x, 0) = U(x, 1) =0 O<x<1, 
UO, y») = U, y)=0 O<y<t, 


where the mesh length is taken to be 4 in each direction. 


D, F 
E D, 


where D,, D are diagonal matrices. Is this ordering consistent with the 
Previous one? How will the asymptotic rates of Convergence for the two 
orderings compare if a solution to the equations is sought using the SOR 
iterative method? 


Given that the eigenvalues of the matrix 


O-o 


are 0, 0, +2, show using the natural ordering of the finite-difference 


a 


` 
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Section B — continued 


Question 10 


Show that for m > 0 and 2 > 0, 


3 1 22 — mè 
[E CIMJ x de = D + r aA 


stating without proof any more basic identities involving Bessel functions 
that you use in the process (J „is the Besscl function of the first kind of order 
m). . 


The equations 


ðu lau, du 
tratz? O0<r<1,0<z<a, 


1,2) =0 O<zca, 
u(r, a) = 0 O<r<l, 
u(r,0) = 1 -r?° O<r<l, 


describe the steady-state temperature u(r, z) in a cylindrical region of 
length a and radius 1 which is insulated along its length and has certain 
radially symmetric temperature distributions maintained at each end. 
Given that the set of functions 


{1, Jo(uar) : Me = the kth positive root of the equation 
Jé) = 0, k = 1,2,...} 


is complete on the interval (0, 1] with respect to the weight function r, 
show that the solution of the problem may formally be expressed as 


a-z < J(u) sinh y,(a — z) 
u(r, 2) = -7 +4 > GT ha oe) - 
kel 
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